Automated Detection of Coronal Loops using a Wavelet Transform Modulus Maxima Method

We propose and test a wavelet transform modulus maxima method for the au- tomated detection and extraction of coronal loops in extreme ultraviolet images of the solar corona. This method decomposes an image into a number of size scales and tracks enhanced power along each ridge corresponding to a coronal loop at each scale. We compare the results across scales and suggest the optimum set of parameters to maximise completeness while minimising detection of noise. For a test coronal image, we compare the global statistics (e.g., number of loops at each length) to previous automated coronal-loop detection algorithms

Solar Magnetic Feature Detection and Tracking for Space Weather Monitoring

We present an automated system for detecting, tracking, and cataloging emerging active regions throughout their evolution and decay using SOHO Michelson Doppler Interferometer (MDI) magnetograms. The SolarMonitor Active Region Tracking (SMART) algorithm relies on consecutive image differencing to remove both quiet-Sun and transient magnetic features, and region-growing techniques to group flux concentrations into classifiable features. We determine magnetic properties such as region size, total flux, flux imbalance, flux emergence rate, Schrijver's R-value, R* (a modified version of R), and Falconer's measurement of non-potentiality. A persistence algorithm is used to associate developed active regions with emerging flux regions in previous measurements, and to track regions beyond the limb through multiple solar rotations. We find that the total number and area of magnetic regions on disk vary with the sunspot cycle. While sunspot numbers are a proxy to the solar magnetic field, SMART offers a direct diagnostic of the surface magnetic field and its variation over timescale of hours to years. SMART will form the basis of the active region extraction and tracking algorithm for the Heliophysics Integrated Observatory (HELIO)

Multiresolution analysis of active region magnetic structure and its correlation with the Mt. Wilson classification and flaring activity

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.Comment: Accepted for publication in Solar Physic

Numerical Simulations of Magnetoacoustic-Gravity Waves in the Solar Atmosphere

We investigate the excitation of magnetoacoustic-gravity waves generated from localized pulses in the gas pressure as well as in vertical component of velocity. These pulses are initially launched at the top of the solar photosphere that is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical and oblique to the gravitational force. We numerically model magnetoacoustic-gravity waves by implementing a realistic (VAL-C) model of solar temperature. We solve two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic-gravity waves, which are later reflected from the transition region due to the large temperature gradient. The atmospheric cavities in the lower solar atmosphere are found to be the ideal places that may act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics

Coronal Mass Ejection Detection using Wavelets, Curvelets and Ridgelets: Applications for Space Weather Monitoring

Coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic feld that can produce adverse space weather at Earth and other locations in the Heliosphere. Due to the intrinsic multiscale nature of features in coronagraph images, wavelet and multiscale image processing techniques are well suited to enhancing the visibility of CMEs and supressing noise. However, wavelets are better suited to identifying point-like features, such as noise or background stars, than to enhancing the visibility of the curved form of a typical CME front. Higher order multiscale techniques, such as ridgelets and curvelets, were therefore explored to characterise the morphology (width, curvature) and kinematics (position, velocity, acceleration) of CMEs. Curvelets in particular were found to be well suited to characterising CME properties in a self-consistent manner. Curvelets are thus likely to be of benefit to autonomous monitoring of CME properties for space weather applications.Comment: Accepted for publication in Advances in Space Research (3 April 2010

An Automated Algorithm to Distinguish and Characterize Solar Flares and Associated Sequential Chromospheric Brightenings

We present a new automated algorithm to identify, track, and characterize small-scale brightening associated with solar eruptive phenomena observed in H{\alpha}. The temporal spatially-localized changes in chromospheric intensities can be separated into two categories: flare ribbons and sequential chromospheric brightenings (SCBs). Within each category of brightening we determine the smallest resolvable locus of pixels, a kernel, and track the temporal evolution of the position and intensity of each kernel. This tracking is accomplished by isolating the eruptive features, identifying kernels, and linking detections between frames into trajectories of kernels. We fully characterize the evolving intensity and morphology of the flare ribbons by observing the tracked flare kernels in aggregate. With the location of SCB and flare kernels identified, they can easily be overlaid on top of complementary data sets to extract Doppler velocities and magnetic field intensities underlying the kernels. This algorithm is adaptable to any dataset to identify and track solar features.Comment: 22 pages, 9 figure

Automated Coronal Hole Detection using Local Intensity Thresholding Techniques

We identify coronal holes using a histogram-based intensity thresholding technique and compare their properties to fast solar wind streams at three different points in the heliosphere. The thresholding technique was tested on EUV and X-ray images obtained using instruments onboard STEREO, SOHO and Hinode. The full-disk images were transformed into Lambert equal-area projection maps and partitioned into a series of overlapping sub-images from which local histograms were extracted. The histograms were used to determine the threshold for the low intensity regions, which were then classified as coronal holes or filaments using magnetograms from the SOHO/MDI. For all three instruments, the local thresholding algorithm was found to successfully determine coronal hole boundaries in a consistent manner. Coronal hole properties extracted using the segmentation algorithm were then compared with in situ measurements of the solar wind at 1 AU from ACE and STEREO. Our results indicate that flux tubes rooted in coronal holes expand super-radially within 1 AU and that larger (smaller) coronal holes result in longer (shorter) duration high-speed solar wind streams

Are Solar Active Regions with Major Flares More Fractal, Multifractal, or Turbulent than Others?

Multiple recent investigations of solar magnetic field measurements have raised claims that the scale-free (fractal) or multiscale (multifractal) parameters inferred from the studied magnetograms may help assess the eruptive potential of solar active regions, or may even help predict major flaring activity stemming from these regions. We investigate these claims here, by testing three widely used scale-free and multiscale parameters, namely, the fractal dimension, the multifractal structure function and its inertial-range exponent, and the turbulent power spectrum and its power-law index, on a comprehensive data set of 370 timeseries of active-region magnetograms (17,733 magnetograms in total) observed by SOHO's Michelson Doppler Imager (MDI) over the entire Solar Cycle 23. We find that both flaring and non-flaring active regions exhibit significant fractality, multifractality, and non-Kolmogorov turbulence but none of the three tested parameters manages to distinguish active regions with major flares from flare-quiet ones. We also find that the multiscale parameters, but not the scale-free fractal dimension, depend sensitively on the spatial resolution and perhaps the observational characteristics of the studied magnetograms. Extending previous works, we attribute the flare-forecasting inability of fractal and multifractal parameters to i) a widespread multiscale complexity caused by a possible underlying self-organization in turbulent solar magnetic structures, flaring and non-flaring alike, and ii) a lack of correlation between the fractal properties of the photosphere and overlying layers, where solar eruptions occur. However useful for understanding solar magnetism, therefore, scale-free and multiscale measures may not be optimal tools for active-region characterization in terms of eruptive ability or, ultimately,for major solar-flare prediction.Comment: 25 pages, 7 figures, 2 tables, Solar Phys., in pres

Multiwavelength studies of MHD waves in the solar chromosphere: An overview of recent results

The chromosphere is a thin layer of the solar atmosphere that bridges the relatively cool photosphere and the intensely heated transition region and corona. Compressible and incompressible waves propagating through the chromosphere can supply significant amounts of energy to the interface region and corona. In recent years an abundance of high-resolution observations from state-of-the-art facilities have provided new and exciting ways of disentangling the characteristics of oscillatory phenomena propagating through the dynamic chromosphere. Coupled with rapid advancements in magnetohydrodynamic wave theory, we are now in an ideal position to thoroughly investigate the role waves play in supplying energy to sustain chromospheric and coronal heating. Here, we review the recent progress made in characterising, categorising and interpreting oscillations manifesting in the solar chromosphere, with an impetus placed on their intrinsic energetics.Comment: 48 pages, 25 figures, accepted into Space Science Review

The Origin, Early Evolution and Predictability of Solar Eruptions

Coronal mass ejections (CMEs) were discovered in the early 1970s when space-borne coronagraphs revealed that eruptions of plasma are ejected from the Sun. Today, it is known that the Sun produces eruptive flares, filament eruptions, coronal mass ejections and failed eruptions; all thought to be due to a release of energy stored in the coronal magnetic field during its drastic reconfiguration. This review discusses the observations and physical mechanisms behind this eruptive activity, with a view to making an assessment of the current capability of forecasting these events for space weather risk and impact mitigation. Whilst a wealth of observations exist, and detailed models have been developed, there still exists a need to draw these approaches together. In particular more realistic models are encouraged in order to asses the full range of complexity of the solar atmosphere and the criteria for which an eruption is formed. From the observational side, a more detailed understanding of the role of photospheric flows and reconnection is needed in order to identify the evolutionary path that ultimately means a magnetic structure will erupt